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Abstract 



Using the CLEO II detector we measure B(Df — ► ne + v) / B(Df — ► (j>e + v) = 
1.24 ± 0.12 ± 0.15, B(D+ n'e+v)/ B(D+ <f>e + v) = 0.43 ± 0.11 ± 0.07 and 
B(D+ r]'e+i/)/B(D+ ije+v) = 0.35 ± 0.09 ± 0.07. We find the vector to 
pseudoscalar ratio, B(D+ 4>e + v)/B(D+ (r) + r)')e + p) = 0.60±0.06±0.06, 
which is similar to the ratio found in non strange D decays. 

PACS numbers: 13.20.Fc, 13.65+i, 14.40.Lb 
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One of the outstanding problems in charm semileptonic decay is the disagreement on 
the vector to pseudoscalar ratio, B(D — > K e + u)/ B(D — > Ke + u). The experimental world 
average of this ratio is 0.56 ±0.06 [1], while theoretical predictions range from 0.5 to 1.2 [1,2]. 
It is of interest to repeat these measurements in the sector, where the initial and final 
hadrons differ by the substitution of a light quark by a strange quark. In the Df sector only 
the Df — > (j)£ + u decay has been extensively studied. The Fermilab experiment E653 has seen 
evidence for the remaining major semileptonic modes, — > (rj + rj')fi + u [3]. In this 
paper we report the first measurements of B(Df — > rje + u) and B(Df — > rj'e + u) normalized 
to B(Df — > (j>e + v), and the vector to pseudoscalar ratio. The data include events with both 
muons and electrons. 

The data consists of an integrated luminosity of 3.11 fb -1 of e + e~ collisions recorded 
with the CLEO II detector at the Cornell Electron Storage Ring (CESR). The data sample 
contains about 3.7 million e + e~ — > cc events taken at center-of-mass energies on the T(45') 
resonance and in the nearby continuum ( y ^/s ~ 10.6 GeV). The CLEO II detector includes a 
Csl electromagnetic calorimeter that provides excellent r\ reconstruction [4]. 

Due to the undetected neutrino, we cannot fully reconstruct Df — > X£ + u decays, where 
X = </>, 7] or 7]' '. However, there are few processes which produce both an X and a lepton 
in the same jet. This correlation is used to extract a clean — > X£ + u signal and is 
implemented by requiring that the X and the lepton be in the same hemisphere with respect 
to the thrust axis of the event. Backgrounds can be reduced for — > rj£ + u and — > (j)£ + v, 
which have sufficient statistics, by also detecting the low energy photon from the D* + — > -D^7 
decay. We refer to this second method as the D* + tag analysis. 

Electron and muon candidates are restricted to lie in the fiducial regions | cos#| < 0.91 
and | cos#| < 0.81, respectively, where 8 is the polar angle of the track with respect to the 
beam axis. Electron candidates are required to have momenta above 0.7 GeV/c and 1.0 
GeV/c for the D* + tag and non D* + tag analyses, respectively. Electrons are identified 
by comparing their ionization energy loss, time-of- flight, energy deposit and shower shapes 
in the electromagnetic calorimeter with that expected for true electrons. Electrons from 
photon conversions and Dalitz decays of 7r°'s are rejected. Muon candidates in the region 
| cos#| < 0.61 (| cos#| > 0.61) are required to have momenta above 1.5 GeV/c (1.9 GeV/c). 
In order to be identified as a muon, a track must penetrate at least 5 interaction lengths of 
iron. 

We identify </>, r\ and r\' candidates by detecting the decays to the K + K~ , 77 and 7/7r + 7r~ 
states, respectively. We require that the momenta of these states be greater than 1.0 GeV/c 
to reduce combinatoric background. Charged kaon (pion) candidates are required to have 
ionization energy loss and time-of- flight consistent with that expected for true kaons (pions). 
Momentum dependent efficiencies for identifying kaons in </> decays are obtained from the data 
by comparing the inclusive </> yield before and after particle identification. These efficiencies 
are then combined with the predicted </> momentum spectrum from — > (j)£ + u Monte Carlo 
events to give the total </> identification efficiency. Photon candidates must lie in the fiducial 
region | cos#| < 0.81 and have a lateral shower shape consistent with that expected for true 
photons. In reconstructing the r\ we require | cos 6d\ < 0.9 to reduce random combinations of 
low momentum photons, where 6d is the photon decay angle in the r\ rest frame with respect 
to the 7] direction in the laboratory. We veto any photon which, when combined with another 
photon, has an invariant mass consistent with the 7r° mass and a momentum greater than 
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0.8 GeV/c. In reconstructing the rj' we require that the invariant mass of the two photons 
from 7] — > 77 decays be consistent with the r\ mass and the r\ momentum be greater than 
0.5 GeV/c. Pion identification efficiency in r\' decays is obtained from the inclusive r\' data 
in the same manner as the kaon identification efficiency for </> reconstruction. 

To suppress the combinatoric background from T(45') events, which tend to be more 
spherical, we require that the ratio of Fox- Wolfram moments [5], _R 2 = H 2 /H , be greater 
than 0.30. The X£ + candidates must have an invariant mass less than 1.9 GeV/c 2 . For 
rj£ + candidates we require that the rj£ + invariant mass be greater than 1.2 GeV/c 2 since 
backgrounds are large at low rj£ + mass. In addition, we require that the X£ + momentum 
be less than 4.5 GeV/c, and greater than 2.0 and 2.5 GeV/c for the D* + tag and non tag 
analyses, respectively. Efficiencies are determined with a Monte Carlo simulation using the 
ISGW model [6]. These Monte Carlo events are passed through a full simulation of the 
CLEO II detector and the same event reconstruction and analysis chain as the data. 

The K + K~ , 77 and r]ir + ir~ invariant mass spectra for all X£ + combinations which pass 
the above selection criteria are fit to obtain the number of candidates which are given in 
Table I. Figures 1(a) and (b) show the invariant mass spectra for 77 and r]ir + ir~. 

There are two main sources of background in the sample: X's accompanied by fake 
leptons [7], and random X£ + combinations. The background due to fake leptons is estimated 
by using the real data to measure the momentum dependent probabilities that a hadron will 
be misidentified as a lepton. These probabilities are typically 0.3% for electrons and 1.2% 
for muons. To estimate the number of events due to fake leptons we take the events which 
pass the basic event selection cuts and treat all charged tracks in these events, which do not 
pass the lepton identification criteria, as leptons. The entries which pass all the cuts are 
weighted by the appropriate fake probability and summed to give the fake background. All 
background estimates and signal yields are given in Table I. 

Random X£ + combinations come from e + e~ — > cc events in which an X produced in the 
fragmentation process is combined with a lepton from a semileptonic decay of the charmed 
hadron in the same jet, and from T(45') — > BB events in which an X and a lepton produced 
in the decay of the B and B mesons are combined. The background from random X£ + 
combinations is estimated using the Monte Carlo simulation. However, the X production 
rate from both fragmentation and B meson decays is not well simulated. For this reason the 
magnitudes of these backgrounds are determined from the data. 

To determine the rate at which fragmentation X's are produced in the same jet as a 
charmed meson we use a sample containing reconstructed D° — > K~ir + decays with D° 
momenta in the same range as those of our sample. Since leptons come primarily from 
charm semileptonic decays, the direction of a charmed hadron is close to that of the high 
momentum lepton. The yields of D° and X mesons produced in the same hemisphere are 
obtained by fitting their invariant mass distributions. In data 1.2 ± 0.5 r\ mesons are found 
for every 1000 reconstructed D° mesons, which is to be compared with 1.7 ± 0.1 in the 
e + e~ — > cc Monte Carlo sample. Therefore for — > rj£ + u we scale the Monte Carlo 
estimate of the charm continuum background by 0.7 ± 0.3. The scale factors for — > (j)£ + u 
and — > 7]'£ + u are 1.0 ± 0.7 and 0.0to!o> respectively. With these factors applied, the 
Monte Carlo simulation predicts the continuum cc background. 

The background from random X£ + combinations in BB events is estimated in a similar 
manner. In this case the directions of the X and the lepton are uncorrelated. Therefore we 



4 



TABLE I. Summary of the non D* + and D* + tag analyses. The errors quoted in this table are 
statistical only. To get the corrected yield, we divide the signal yield by the detection efficiency (e) 
and the branching fraction (£>) for the subsequent decays of <f>, tj or tj' mesons, and a factor of two 
to take into account the fact that the lepton can be either electron or muon. 





non D* + tag 


D *+ 


tag 


Decay mode 


D+ <j>l+v 


D+ T)£ + P 


D+ rj'l+v 


D+ <j>l + v 


D+ r]l + v 


Candidates 


863.6 ±40.5 


577.6 ±30.3 


42.3 ± 7.9 


326.9 ±27.1 


153.4 ± 15.0 


Fake lepton 


95.3 ±2.2 


128.8 ± 1.6 


9.0 ± 0.4 


27.9 ± 1.2 


22.7 ± 0.7 


cc -\- JD JD 


72.7 ± 2.0 


62.7 ± 3.2 


1.9 ± 0.2 


Zo.O ZL l.U 


1 o Q J_ 1 A 
iz.o m i.^t 


D+ -> 7]'£ + V 




5.0 ± 1.3 






2.3 ± 0.6 


D+ + D+ 


695.6 ±40.6 


381.0 ±30.5 


31.4 ± 7.9 






fc 




0.855 


0.926 




0.989 


Df signal yield 


695.6 ±40.6 


325.6 ±26.1 


29.1 ± 7.3 






D*+ + Photon bkgd. 








275.5 ±27.2 


114.8 ± 15.1 


fsig 








0.759 


0.803 


D* + signal yield 








209.2 ±20.6 


92.2 ± 12.1 


e • B(%) 


2.078 


0.807 


0.204 


1.408 


0.474 


Corrected yield 


16734 ± 976 


20184 ± 1616 


7128 ± 1786 


7431 ± 733 


9730 ± 1276 



compare the number of X's with momentum above 1.0 GeV/c in the continuum subtracted 
T(45') data with the number observed in the T(45') BB Monte Carlo sample. The resulting 
backgrounds are found by scaling the Monte Carlo estimate. 

For Df — > rj£ + u decay, we estimate the Df — > rj'£ + u feed down by multiplying the 
efficiency corrected yield of the Df — > rj'£ + u decay in our data sample, as discussed below, by 
the appropriate efficiency and the branching fraction. The Cabibbo suppressed D + — > rj£ + u 
decay contributes a sizable contamination to the Df — > rj£ + u decay in the non D* + tag 
analysis. Using a measurement [13] of the ratio of a(D + )/ a(Df) and the ISGW2 model 
[10] for the ratio of widths, T(D + — > rj(r] l )£ + u)/T(Df — > r]{ri')£ + v), we estimate the fraction 
of the yield which is from D+ to be = 0.855 ± 0.051(0.926 ± 0.026) for an rj - rj' 

mixing angle of ( — 15 ± 5)°. 

Figures 1(c) and (d) show the fitted number of 7/'s and rj n s which fall in each bin of the 
r}£ + and rj'£ + mass, respectively. The combined background estimate is also shown, as well 
as the sum of the simulated signal and the combined background estimate. The simulated 
signal has been normalized to the number of candidates extracted from the fit to the 77 and 
7/7r + 7r~ invariant mass spectra. Table I gives the yields after correcting for efficiencies and 
the (j) — > K + K~ , 7] — > 77 and r\' — > r]ir + ir~ branching fractions. 

In the D* + tag analysis, we look for a low energy photon from the D* + — > Df^ decay 
in the same jet as the X(= (j) or rj) and the lepton. The photon candidate must now 
lie in the fiducial region | cos#| < 0.71, and have energy greater than 0.12 GeV. To further 
suppress background photons from 7r° decays, we veto any photon that has an invariant mass 
consistent with the 7r° mass when combined with another photon. To select X£ + ~f candidates 



5 



we require that AM = Mj<+ 7 — Mxi+ be between 0.1 and 0.2 GeV/c 2 , where Mxi+ and 
Mxi+j are the invariant masses of X£ + and X£ + ~f systems, respectively. Figure 2(a) and (b) 
show the K + K~ and 77 invariant mass spectra for all (j>£ + ^ and rj£ + ^ combinations which 
pass the above selection criteria. 

The backgrounds due to fake leptons, random X£ + combinations, and the Df — > rj'£ + u 
feed down are estimated by the same procedures as for the non D* + tag analysis. The 
contamination from the Cabibbo suppressed decay D + — > rj£ + u is estimated to be small. The 
background from the combination of a true X£ + and a random photon is estimated using 
a Monte Carlo simulation. To test the agreement between the Monte Carlo estimate and 
the data for this background, we use events with the decay chain, D* + — > D°ir + and D° — > 
K~e + u. These events are ideally suited for the test because a K~e + pair has kinematics 
similar to that of an X£ + pair from the decay. We combine photons with K~e + pairs 
selected with the same criteria as for X£ + 's from the decay. In data 21.4 ± 1.7 if~e + 7's 
are found for every 100 reconstructed K~e + pairs, which is to be compared with 22.3±0.1 for 
the Monte Carlo sample. Therefore we scale the Monte Carlo estimate of the random photon 
background by the factor of 0.96 ± 0.08. With this scale factor, the Monte Carlo simulation 
predicts the ratio of the signal to the sum of the signal and the random photon background, 
f sig , to be 0.759 ± 0.030 for D+ -> <f>£+u and 0.803 ± 0.031 for D+ -> rj£+u. Figures 2(c) 
and (d) show the fitted number of </>'s and 7/'s which fall in each AM bin, respectively. The 
combined background estimate is also shown, as well as the sum of the simulated signal 
and the combined background estimate. The simulated signal has been normalized to the 
number of candidates extracted from the fit to the K + K~ and 77 invariant mass spectra. 

To obtain the effective branching ratio in the electron channel, we take into account the 
reduced phase space for the muon channel [8], and the efficiency loss due to the final state 
radiation for the electron channel [9]. The results for B(Df — > rje + u)/ B(Df — > (j>e + u) are 
1.21 ± 0.12 ± 0.16 and 1.32 ± 0.22 ± 0.15 for the nontag and tag analyses, respectively. The 
two measurements for r)/<f> are combined taking into account the correlated statistical and 
systematic errors. This result and the r)' / <f> measurement from the nontag analysis are given 
in Table II. 

The major contributions to the systematic error in B(Df — > rje + u)/ B(Df — > (j>e + u) 
are the uncertainties in the D + — > rj£ + u contamination, fitting, fake lepton, and charm 
continuum backgrounds. The systematic errors in B(Df — > rj'e + u)/B(Df — > (j>e + u) and 
B(Df — > rj'e + u)/B(Df — > rje + u) are dominated by the uncertainty in the charm continuum 
background. 

Model predictions and measurements are listed in Table II. Our measurements of 
B(Df -> rje+u)/B(D+ -> <f>e + u) and B(Df -> rj'e + u)/B(D+ -> <f>e+u) agree well with 
the ISGW2 model predictions with an r\ — r\' mixing angle of —10° [10]. Our measure- 
ment of the vector to pseudoscalar ratio for the Df semileptonic decays, 0.60 ± 0.06 ± 0.06, 
agrees well with the measurements for D semileptonic decays, 0.56 ± 0.06 [1], and also with 
the ISGW2 model prediction [10]. Using the factorization hypothesis, Kamal et al. pre- 
dict B(D+ -»■ 7]'e+u)/B(Dt -> rje+u) = B(D+ -»■ r]'p + )IB{Dt ->■ VP + ) [H]- We obtain 
B(D+ -> ri'p + )/B(Df -> V p+) = 1.20 ± 0.35 from the CLEO measurements [12]. This dis- 
agrees with our measurement of B(Df — > rj'e + u)/ B(Df — > rje + u) = 0.35 ± 0.09 ± 0.07 from 
the nontag analysis. 

In conclusion, we have measured B(Df — > rje + u)/ B(Df — > (j>e + u) = 1.24 ± 0.12 ± 0.15 
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TABLE II. Summary of measurements and predictions. The numbers in (parentheses) [brack- 
ets] are the ISGW2 predictions using values of the 77 — 77' mixing angle of ( — 10°) and [—20°]. 



This Experiment 



E653 



ISGW2 



B(D+^4>e+v) 
B(Dj^4>e+v) 



0.43 ±0.11 ±0.07 



0.60 ±0.06 ±0.06 



1.24 ±0.12 ±0.15 



< 1.6 @90%C.L. 




(0.50)[0.67] 



(0.60)[0.69] 



(1.17)[0.77] 



B{D+^{r 1 +r 1 ')e+ V ) 



and B(Df -> rj'e + u)/B(Df -> <f>e+u) = 0.43 ± 0.11 ± 0.07. Our value for the ratio of vector 
to pseudoscalar rates for the Df semileptonic decays agrees with that observed for the D 
semileptonic decays, and also with the ISGW2 model prediction. This agreement increases 
confidence in the evaluations of the Df — > (j)ir + branching ratio using models of semileptonic 
decays and measurements of T(Df — > (j)£ + u)/T(Df — > (j)ir + ) [1,2]. 
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luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Heisenberg Foundation, the Alexander von 
Humboldt Stiftung, the Natural Sciences and Engineering Research Council of Canada, and 
the A. P. Sloan Foundation. 
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FIG. 1. The (a) 77 and (b) tjtt + tt~ invariant mass spectra for _D+ — ► tj£ + i/ and _D+ — ► tj'£ + i/ 
candidates in the non D* + tag analysis: The solid curve is a fit to each spectrum. The (c) tj£ + 
and (d) tj'£ + invariant mass spectra for the candidates: The points with error bars represent the 
number of candidates in each mass bin. The solid histogram shows the simulated signal plus the 
predicted background, and the dotted histogram shows the predicted background. We removed the 
cut on T}£ + and tj'£ + invariant masses discussed in the text for (c) and (d). 



FIG. 2. The (a) K + K~ and (b) 77 invariant mass spectra for Df — ► (f>l + v and Df — ► tj£ + i/ 
candidates in the D* + tag analysis: Note that a candidate photon from the D* + is required. 
The solid curve is a fit to each spectrum. The spectra of the pseudo mass difference AM for (c) 
_D+ — ► (j)£ + i/ and (d) _D+ — ► tj£ + i/ candidates: The points with error bars represent the number 
of candidates in each AM bin. The solid histogram shows the simulated signal plus the predicted 
background, and the dotted histogram shows the predicted background. We removed the cut on 
AM discussed in the text for (c) and (d). 
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